On perfect fluid models in non-comoving observational spherical coordinates 
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We use null spherical (observational) coordinates to describe a class of inhomogeneous cosmolog- 
ical models. The proposed cosmological construction is based on the observer past null cone. A 
known difficulty in using inhomogeneous models is that the null geodesic equation is not integrable 
in general. Our choice of null coordinates solves the radial ingoing null geodesic by construction. 
Furthermore, we use an approach where the velocity field is uniquely calculated from the metric 
rather than put in by hand. Conveniently, this allows us to explore models in a non-comoving frame 
of reference. In this frame, we find that the velocity field has shear, acceleration and expansion rate 
in general. We show that a comoving frame is not compatible with expanding perfect fluid mod- 
els in the coordinates proposed and dust models are simply not possible. We describe the models 
in a non-comoving frame. We use the dust models in a non-comoving frame to outline a fitting 
procedure. 

PACS numbers: 04.20.Cv, 04.20.Jb, 98.80.Jk, 04.40.Dg 



I. INTRODUCTION 

The study of inhomogeneous cosmological models is 
a well motivated and justified endeavor (see for reviews 
^2]). These models provide more freedom in discussing 
very early or very late evolution of the irregularities in the 
Universe. Their study also complements perturbation 
approaches. It is worth mentioning that there are a few 
hundreds of inhomogeneous cosmological models that re- 
produce a metric of the Friedmann-Lemaitre- Robertson- 
Walker (FLRW) class of solutions when their arbitrary 
constants or functions take certain limiting values 
They become then, in that limit, compatible with the 
almost homogeneous and almost isotropic observed Uni- 
verse. This shows the richness of these studies. 

A difficulty that is encountered in this models is that 
the null geodesic equation is not integrable in general. In 
this paper, we explore the alternative of using null (ob- 
servational) spherical coordinates in which the radial null 
geodesic equation of interest is solved by construction. 
However, when considering null coordinates and a given 
metric for the spacetime some subtleties arise regarding 
the frame of reference used. In order to explore this point 
we will use in this paper the approach described in Ishak 
and Lake0| where the velocity field is calculated from 
the metric and not put in by hand. Conveniently, this 
approach allows one to explore non-comoving frames of 
reference, an important point for this paper. 

Surprisingly, little work has been done in non- 
comoving coordinates 0, H 0, H H in im despite some 
interesting features particular to them. Notably, there 
are models which are separable only in a non-comoving 
coordinate system p^j. Moreover, exact solutions to Ein- 
stein's equations in a non-comoving frame usually have a 
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rich kinematics with shear, acceleration and expansion. 
Such solutions are relatively rare in the comoving frame 
Isl , see also a recent discussion in Till . Another point dis- 
cussed in [To| is that comoving coordinates do not cover 
all the spacetime manifold for a specified energy momen- 
tum tensor. Finally, it is worth mentioning that it is 
often difficult to do the mathematical transformation of 
a given solution from non-comoving coordinates to co- 
moving ones, and even when the passage is made, there 
is no guarantee that the solution will continue to have a 
simple or explicit form. 

In the present paper, we use the spherical null Bondi 
metric | i25j to present models where a non-comoving 
frame is proven necessary. We explicitly demonstrate 
how a comoving frame leads to severe limitations. 

Furthermore, we use the dust models in the non- 
comoving frame to outline a fitting procedure where ob- 
servational data can be used to integrate explicitly for 
the metric functions. Using observational coordinates is 
particularly useful when one wants to compare directly 
an inhomogeneous model to observational data. Such an 
interesting progra m had been nicely developed in Refs. 
m, m 13, [IE IS El where the authors used a general 
metric that can be written as a FLRW metric plus exact 
perturbations. The spherically symmetric dust solutions 
were considered in Ref . ■ The authors also developed 
and used a fluid-ray tetrad formalism 13] in order to de- 
rive a fitting procedure where observations can be used 
to solve the Einstein's field equations. After some neces- 
sary revisions , this program has been relaunched 
recently [Ilial"^ 

We consider here in our work the spherically symmetric 
case but using the Bondi metric 25] in a non-comoving 
frame. Also, we don't use the fluid-ray tetrad formal- 
ism |13| but the inverse approach to Einstein's equations 
developed in 0. 

In the following section, we set the notation and recall 
some useful results. In section III, we discuss observa- 
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tional coordinates and explain the cosmological construc- 
tion around our world-line. We also discuss the physical 
meaning of the functions that appear in the metric used 
here. We provide in section IV perfect fluid models in 
a non-comoving frame. In section V, we show how dust 
models are not possible in a comoving frame. We describe 
dust models in a non-comoving frame and outline a fit- 
ting procedure in section VI and summarize in Section 
VII. 



Pi and p2 are the pressures respectively parallel and per- 
pendicular to n". When pi = p2 and the shear term 
vanishes the fluid is called perfect. It was shown in 
that in the case where 



A = a^n°'nf^ ^ 0, 



we have 



Gl 



277 ' 



(6) 
(7) 



II. NOTATION AND PRELIMINARIES 

We set here the notation and summarize results to be 
used in this paper. In Ref. warped product space- 
times of class Bi [21I l2^ were considered. These can be 
written in the form 



(1) 



where C(x") ~ r{x^ , x^yw{x^ , x^)"^ , sig{Tii) = and 
s«(/(E2) ^ 2e {e ~ ±1). Although very special, these 
spaces include many of interest, for example, all spheri- 
cal, plane, and hyperbolic spacetimes. For Si, we write 



a{dx^f + 2bdx^dx^ + c{dx'^f 



(2) 



with a,b and c functions of only. Consider a 

congruence of unit timelike vectors (velocity field) = 
(m^,m^,0, 0) with an associated unit normal field (in 
the tangent space of Ei) satisfying Uau" = 0,nan°' = 1 
|23j . It was shown in 4] that u" is uniquely determined 
from the zero flux condition 



(3) 



where is the Einstein tensor of the spacetime. The 
explicit forms for and were written out for canonical 
representations of Ei, including the null (Bondi) type 
of coordinates that we use in the present paper. With 
G = G^, Gl = G^w"u0 and G2 = G^n^np, it was 
shown in H that the condition 



G -I- Gl = 3G2 



(4) 



is a necessary condition for a perfect fluid source, and 
that in some cases, this condition is also sufficient. For 
example, in |^, equation ^ was used to derive an algo- 
rithm which generates all regular static spherically sym- 
metric perfect fluid solutions of Einstein's equations. In 
this paper, we are interested in perfect fluid sources so 
it is important to recapitulate the following results from 
. Consider a fluid with anisotropic pressure and shear 
viscosity but zero energy flux (non-conducting). The 
energy-momentum reads 

Tp = pu"ui3+pin°'np+p2Sp+p2{u°'ui3-n"nf3)-2ria'0, 

(5) 

where p is the energy density and is the shear associ- 
ated with u"; ?7 is the phenomenological shear viscosity; 



pi 



G2 
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G + Gl - G2 





(8) 



(9) 



and 77 is a freely specified function. The procedure to 
impose a perfect fiuid source in this degenerate case is to 
impose the condition Q and also necessarily set ?/ = 0. 
For other choices of 77, the fluid is viscous. 



III. THE METRIC AND OBSERVATIONAL 
COORDINATES 

We consider in the present paper the null coordinate 
system {x""} — {v, r, 9, </)}. These are called observational 
(or cosmological) coordinates as we can construct them 
around our galaxy world line C as indicated in Figure 
n The trajectory defined by v, 9 and (p constant is a 
radial null geodesic and each hyper-surface of constant 
z; is a past light cone of events on C. We choose v — Vo 
and r = to represent the vertex "here and now" . The 
coordinate r is then set by construction to be the area 
distance as explained further and is related to the lu- 
minosity distance by r = d^/{l + z)^ (see e.g.0). 
Finally, 9 and (p are the spherical coordinates on the ce- 
lestial sphere. The geometry of the models is represented 
by the general spherical Bondi metric in advanced coor- 
dinates t2S |: 



ds ~ 2cdr dv 



o , 2m , T 

1 )dv^ 

r 



r^{d9'^ + sin{9fd(j)'^) 
(10) 

where c = c(r, u) > 0, m = m(r, v) and r > 2m. The 
radial (9 and constant) ingoing null geodesic equation 
w=constant is solved by construction (see Appendix 0. 
The components of the mixed Einstein tensor G^ for H1(J|) 
are given in appendix ^ and the structure of the Weyl 
tensor is discussed in Appendix |0 Regularity of the 
metric and the Weyl invariants requires that m(r, v) and 
c(r, v) are G^ at r = (e.g. see equation ijnjj). It follows 
that 



(1 



2m(r, v) 



1. 



(11) 



Also, we can use the freedom in the null coordinate v 
to normalize it by setting c(0,u) = 1. As we will write 
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witmvo, r,theta2,[)hi2) 



is the null 4-vector tangent to the congruence of null 
geodesies. The physical meaning of the 0£tical scalars 
can be understood in the following way [3^ l3^ : If an 
opaque object is displaced an infinitesimal distance dr 
from a screen (perpendicularly to the beam of light), 
it will cast on the screen a shadow that is expanded 
by Oopticaidr, rotated by ujopUcaidr and sheared by 
\(Topticai\dr. As expected from the spherical symmetry 
of the geometry, the non-vanishing optical scalar for the 
null congruence is the optical rate of expansion, from 
which wc find 

rc{r,v)^—^ (17) 

"optical 

We identify from equation H17|) that rc(r, v) is a measure 
of the reciprocal of the expansion of null rays. 



IV. MODELS IN A NON-COMOVING FRAME 



FIG. 1: Observational coordinates {v,r,0,(l)}. Our past null 
cone is defined hy v = Vo- C is the observer world-line, C 
is the world-line of another celestial object. The trajectory 
defined by v, 6 and (j) constant is a radial null geodesic, u" is 
the fluid velocity field vector and k" is the null tangent vector. 
The null rays are traveling in the opposite direction to k" as 
drawn on the flgure. The coordinate r increases along the 
trajectory CC down the light-cone. 



further in this paper (see equation (|52ll ') this means that 
we require that v measures the proper time r along our 
galaxy world line C. 

Whereas the meaning of the metric function m(r, v) 
is very well known, we are not aware of any previous 
literature where an interpretation for c(r, v) was given. 
The function m(r, v) represents the effective gravitational 
(geometrical) mass (e.g. QJ\, ^J, _29J, ^ and [sj) and 
is given by 



3/2 



-Ra 



(12) 



where Re^,^'^ is the mixed angular component of the Rie- 
mann curvature tensor. For the physical meaning of the 
function c{r,v), it turns out to be useful to study the 
kinematics of null rays. These usually include the op- 
tical shear, vorticity and rate of expansion, respectively 
defined by [13 



A. The velocity field 

We consider an observer moving with a fluid for which 
the streamlines are given by the general radial timelike 
vector u" — [u^{r, v),u'^{r, v), 0, 0]. We assume that such 
a velocity field exists for which the energy-momentum 
tensor takes the perfect fluid form 



(18) 



where p and p are respectively the energy-density and 
isotropic pressure associated with u". The velocity field 
is simply determined from the zero fiux condition (jSJ and 
is given by 



1 



1 



(19) 



u — u 



where ' = and * 

or 



C{r,v) \ (l _ 2m{r,v) \2 4m'(r,v) 

V V r J rc' (r.v) 

2m{r,v)\ n { A 2m(r,u)\2 4m*(r,t;)\ 

V~ r )~ r j + Tc'(r.v) I 

^ (\ 2m(r,i!)\2 4m»(r,v) 

Y V r } rc' {r,v) 

(20) 

The associated unit normal 



dv ' 



vector field n" (uau" — and Uan" = 1) is given by 

n„ = c(r,'^;)K,-u^O,0]. (21) 



2 = 
optical — 

optical 




(13) 
(14) 


Interestingly, the velocity field has shear ct^ 
celeration 7^ and expansion rate scalar ( 
general. 


7^ 0, ac- 
9 7^ in 


optical = 


2 


(15) 


B. The perfect fluid condition 





where 



fc„ = [1,0,0,0] 



It follows from the metric (|10|l that A, as defined in 
(16) ©, is not zero and that for a perfect fluid source we 
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must impose the condition Q and set 77 = 0. With 
m = m(r, v) and c = c(r, w), condition reads 



£ - c^V^ = 



(22) 



where 



and 



C = c^{2m' -m"r) + c'^{3c'{m-rm') (23) 
+c"r(r - 2m)) + cr^c" - cV^c* 



Af = c'{r ~ 2vif + 4rm*. 



(24) 



The metric (jlOf) along with the metric constraint H22f) 
represent a perfect fluid model with 



Gl 2(cm)'-c'r + 



and 



P (=Pl =P2) 



-2{cm)' + cV 



(25) 



(26) 



V. MODELS IN A COMOVING FRAME 

In this section, we specialize to models in a comoving 
frame of reference. We show how this frame fails in the 
realization of the cosmological construction proposed. 



for a perfect fluid model is equation Q), which can be 
written as 

- c^cV + bc^c'm + 2c^m' + c^c'V^ - 2c^c"mr- 

Sc^c'm'r - c^m"r - c'c'r'^ + d'cr"^ = 0. (32) 

For a perfect fluid source in this frame the pressure p is a 
function of both r and v while the energy density is only 
a function of r. 



Gl _ m'{r) 
Stt 47rr2 



(33) 



and 



c'(r, v){r + 2m(r)) + c(r, v)m'(r) ,„ ^, 

P = A — TT ^ ■ (^4'' 

47rr^c(r, v) 

The 4-acceleration it" has the non- vanishing components 
. ^ rc' [r, v)(r — 2m{r)) + c(m(r) — m! (r)r) 



rc'^{r,v){r — 2m(r)) 



(35) 



and 



rc'{r,v){r — 2m{r)) + c{m{r) — m'{r)r) 
u — 5—^ r . (3d) 

r^c(r, w) 

A caveat in this frame (comoving) is that the expansion 
scalar vanishes and the model is not suitable for describ- 
ing an expanding Universe. 



A. Perfect fluid models 



B. The zero-pressure case 



With the metric function 



-c\r,v)(l 



2m(r, v) 



) < 0, (27) 



the requirement of comoving coordinates reads 

dm{r,v) 

u"" = ^ Gi: = 2 = 0. 



(28) 



Hence, the necessary and sufficient condition for a co- 
moving frame is 



m(r, v) = m(r). 
It follows from and that 



c(r,.)^(l-H!IlM) 



1- 



2m{r) 



(29) 



(30) 



(31) 



and Uy = 0. With this velocity field the shear tensor 
vanishes, therefore, the necessary and sufficient condition 



The present matter dominated (as opposed to radia- 
tion dominated) Universe is well approximated by a zero- 
pressure model, commonly referred to as "dust" . In this 
case (comoving) A = and the zero pressure conditions 
follow from equations (|SJ) and (O as 



G2 = 



and 



G + Gl = 0. 
With (EHJ, equations and (jSHJ read 



c r — 2m 



= 0, 



(37) 
(38) 

(39) 



and 



2 / , 2 f ,2/12 n, 2 II 2 I I 

ccr + ccm + ccr — Zc c rm — 6c c rm — 

c m r — c c r + c cr = v. (4U) 



Integrating gives 



^-ir,v) = /(w)exp 



m'{r) 
r — 2m(r) 



dr) (41) 
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which when put into H40() gives 



r^{r — 2m(r)) 



0. 



(42) 



With f{v) > (from c(r, v) > 0) the zero pressure model 
reduces to the following two cases: 

i) m(r) = and the spacetime reduces to the Minkowski 
flat spacetime {RaPfS = 0), or, 

ii) m'(r) = (m is constant) and the spacetime reduces 
to Schwarzschild vacuum in Eddington-Finkelstein coor- 
dinates {RajS = 0, Rap-yS 7^ 0.) 

Therefore, a dust model is not possible in a comoving 
frame using the observational coordinates and the Bondi 
metric We turn in the following section to a non- 

comoving frame for dust models. 



VI. DUST MODELS IN A NON-COMOVING 
FRAME 

A. The velocity field 

We are interested in building dust models using spheri- 
cal observational coordinates and a non-comoving frame. 
In a 1+3 decomposition of the spacetime, these models 
are given by the general Lemaitre-Tolman-Bondi solution 
[33 | [j] . In this non-comoving case A ^ in general, so we 
must set ry = and impose the zero pressure conditions 
(|37|l and H38() which can be written as 



,/CTO / 

TO = cm 

V cr V 

and H38|l can be written as 



- 1- 



2to, 



;(1(3cW + cto")-^(1 + ^)-c"(1-^; 



(43) 



c c 
(44) 



The metric (|10|) with constraints (|43|l and (|44|l represents 
a class of inhomogeneous dust models in spherical ob- 
servational coordinates. Using equation 0, the energy 
density is given by 



At: p{r, v) 



2m'{r,v) c'{r,v) / 2m{r,v) 



rc{r, v) 



,«) V 



(45) 



This result can also be obtained from the effective grav- 
itational mass equation (|12|l . The velocity field follows 
from equations (|19|l and H2Q(I : 



1 



c(r, v) 



(1- 



2m{r.v) 



) + 



2m* {r.v) 
m' c 



(46) 



or equivalently by using H43|l 
1 



c(r, v) 



2m' (r,v)c(r,v) 
rc' (r,v) 



(1 



2m(r,v) 



(47) 



and 



m*{r, v) 


m'(r, v) 



(48) 



We verified that the acceleration 4-vector field it" van- 
ishes as the dust fluid is moving geodesically. Interest- 
ingly, the velocity field remains with non- vanishing shear 
and expansion rate. 



B. The conformally flat case 

It is a well known result that a cosmological model 
which satisfies the Einstein equations with a perfect fluid 
source, a barotropic equation of state i.e. p — p{p) (in- 
cluding p = 0), which is conformally fiat {Cafs-yS = 0) 
and has non-zero expansion is a Lemaitre-Friedmann- 
Robertson- Walker model (LFRW) Q, and For 
dust models in the non-comoving frame, the condition 
Ca/37(5 = (see appendix Q reduces to 



c' / 2m \ 'Zm' 'dm 

-(1 + — 

c \ r / r T 



2m' 
r 



3to 



= 



(49) 



Therefore, the metric HlOjl along with the constraints 
(|43|l . H44|) and H49I) represents the homogeneous and 
isotropic (LFRW) limit of the models. We are interested 
here in more general inhomogeneous models. 



C. Basic observable quantities 

1. The Redshift 

The light emitted with a wavelength Ae from a point on 
the light cone is observed at the vertex "here and now" 
(see Figure ^) with a wavelength Aq. The redshift is then 
given by (see e.g. [s^, 



l + z 



e 



'hserver 



{kpuf^)obse 



A. 



(50) 



where u" is the normalized timelike velocity vector field 
and ka is the null vector as given previously by Hlt)|l . It 
follows that 



kau" = u"{r, Vo) 



(51) 



where u'"{r, Vo) is given by l|19|l . It follows from the regu- 
larity condition Hll|l and the timelike normalization con- 
dition m"mq = — 1 evaluated at r = (observer) that 



= zi''(0,^;o) = — ^ 

observer ^y^-) '^o J 



1 



(52) 



where in the last step we used the freedom in the null 
coordinate v to set c(0, v) = 0. Finally, 



l + z 



^ '^o^ emitter 
'ii^(r, Vo) observer 



u'^{r, Vo) emitter • 

(53) 
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For the dust case, equation gives 
1 1 



(54) 



or equivalently by using the constraint H43|) 



c(r, Vq) / m'{r,Vo)c(r Vq) 7n*{r,Vo) ~ 

'* rc' {r,V(y) m' {r,Vo)c{r,Vo) 



(55) 



where we have set r = remitter in equations H54(l and (|55|l . 



S. The observer area distance 

The coordinate r in the model is set by construction 
to be the observer area distance [s^ which is defined 
by dA = r^dVl where dA is the cross-sectional area of 
an emitting object, and dO. is the solid angle subtended 
by that object at the observer. The area distance r is 
related to the luminosity distance dL hy r = d^/ (\ + 
zY [i^' The luminosity distance can be determined by 
comparing the observed luminosity of an object to its 
known intrinsic luminosity: 47rd^ = L/F where F is the 
observed (measured) flux of light received and L is the 
object's intrinsic luminosity. 



3. Galaxy number counts 

Another observable of interest is the source number 
counts as a function of the redshift. An observer at the 
vertex "here and now" will count on the light cone a 
number dN of sources between redshifts z and z + rfz in 
a solid angle dU.. It follows that 



—r- = Jcn(Vo,r)r dil — 
dz dz 



(56) 



where n{vo, r) is the number density of sources and fc is a 
fractional number indicating the efficiency of the counts 
(completeness) |0. This number corrects for errors in 
source selection and detection, see e.g. 0, |3^. For 
simplicity, we can assume that the necessary adjustments 
for the dark matter can be incorporated via fc- The 
energy-density follows 



p{vo,r) = n{vo,r)M 



(57) 



where M is the average rest mass for the counted sources. 



D. A fitting procedure algorithm 

As discussed earlier, the approach used allows us to 
integrate explicitly the model given observational data. 
As a first step, we rearrange the model equations. We 



combine equations H45() evaluated at v = Vo with equation 
and use c(0, Vq) = 1 to obtain 



C'[r, Vq) ^ 

which integrates to 
c(r, Vo) 



A7T{l + z)^{r,Vo)pir,Vo)) (58) 



(59) 



Y^l - 8 TT /(I + z)^ (r, Vo)rp{r, Vo)dr 



Integrating H45|) for ■m{r,Vo) gives 
1 



m(r, Vo) 



2c(r, Vo) 



c'{r, Vo) + 47rrp(r, Vo)c{r, Vo)j rdr 

(60) 

where we also used m{0,Vo) — 0. Now, the observations 
provided as polynomial functions p{z) and r(z) fitted to 
the data can be used to integrate explicitly for m{r,v) 
and c{r,v). The steps for the fitting algorithm are as 
follows: 

• Express cosmological data as polynomial functions 
for: 

i) The energy-density p{r, Vo) from galaxy number 
counts. Many projects are accumulating very large 
amounts of data. See, for example, [39j for the 
"Sloan Digital Sky Survey" . 

ii) The observer area distance r(z, Vo) from "stan- 
dard candles" projects in which it is possible 
to measure the redshift and the distance in- 
dependently. The accumulating data from the 
supernovae cosmology projects are very promis- 
ing. See, for example, [43 for the "High-Z SN 
Search" project, 41] for the "Supernova Cosmology 
Project" and 42J for "The Supernova Acceleration 
Probe (SNAP)" project. 

• Invert the function r{z,Vo) to obtain z{r,Vo)- 

• This can in turn be used to write the energy-density 
polynomial function as p{z{r, Vo), Vo)- 

• Now, with z(r, Vo) and p(r, Vo) expressed as func- 
tions of r (and not z), integrate equation (|59|l over 
r to obtain c(r, Vo) on the light cone. 

• With c(r,Vo) determined, integrate equation H6U|) 
over r to obtain TO(r, Vo) on the light cone. 

• Finally, with c{r,Vo) and m{r,Vo) determined, use 
equations H43|) and (|44|l to integrate over v. 

The level of difficulty of this last step can be monitored 
using the analytical forms used for p(z) and r(z) and 
remains a tractable problem while integrating the null 
geodesic equation in the standard 1-1-3 form of the LTB 
models is not tractable (see e.g. llSI') . and one has to 
recourse to numerical integrations |43. Moreover, the 
fitting procedure has the interesting feature of incorpo- 
rating the observations in the process of integrating ex- 
plicitly for the metric functions. 
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It is worth mentioning that in principle the information 
on our hght cone can not determine its future evolution 
uniquely. We need to make the reasonable assumption 
that there will be no future events in the cosmic evolu- 
tion that will invalidate the entire data obtained from our 
light cone, see e.g. Furthermore, one must keep in 

mind the usual limitation of the underlying models used 
here as they are spherically symmetric around our world- 
line and more general inhomogeneous models should be 
considered in future studies of fitting procedures. 



For the Bondi metric H10() . the 4 equations I|A1|I are all 
satisfied for v = constant. 

APPENDIX B: EXPRESSIONS FOR Gp 
COMPONENTS 

The expressions for the components of the mixed Ein- 
stein tensor are as follows 

cr J. ' j.^ 



VII. SUMMARY 



GV 



2c' {r, v) 



(B2) 



We expressed inhomogeneous cosmological models in 
null spherical non-comoving coordinates using the Bondi 
spherical metric. A known difficulty in using inhomoge- 
neous models is that the null geodesic equation is not 
integrable in general. Our choice of null coordinates 
solves the radial null geodesic by construction. We iden- 
tified the general meaning of the metric function c(r, v) 
to be the reciprocal of the optical expansion. We used 
an approach where the velocity field is uniquely calcu- 
lated from the metric rather than put in by hand. Con- 
veniently, this allowed us to explore models in a non- 
comoving frame of reference. In this frame, we find that 
the velocity field has shear, acceleration and expansion 
rate in general. In this set of coordinates, we showed 
that a comoving frame is not compatible with expanding 
perfect fluid models and dust models are simply not pos- 
sible in this frame. We described then perfect fluid and 
dust models in a non-comoving frame. The framework 
developed allows one to outline a fitting procedure where 
observational data can be used directly to integrate ex- 
plicitly for the models. 
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APPENDIX A: NULL GEODESIC EQUATION 

The paths of light rays are described by null geodesic 
trajectories under the eikonal assumption 26|. The 
geodesic trajectories are determined by solving the null 
geodesic equation 



«q;/3 



(Al) 



where fc" is a null vector {k^ka — 0) tangent to the 



Gr 



GV 



2m* 



2m' 

^2 



cr 
d 



2m \ m" c" 
— + — 



(B3) 

(B4) 

(B5) 
(B6) 



where * = ' = c = c(r, w) and m = m{r,v). We 
note that these components are related by 



(B7) 



APPENDIX C: THE WEYL TENSOR AND THE 
CONDITION FOR CONFORMAL FLATNESS 

The structure of the Weyl tensor Caf3jS = is usually 
explored to derive the conformally flat case of a cosmo- 
logical solution (i.e. Cap-yS = 0). This can reveal the 
limits of the model's parameters for which it reduces to a 
Lemaitre-Friedmann-Robertson- Walker (LFRW) model. 
The non- vanishing components of the mixed Weyl tensor 
for the metric (|10|) are related and given by 



(CI) 



where 



2m 

r 



2a/'' = 2a/*->V(r,t;) 

cV(m"r - 4m' + —) - c^^'r"" (C2) 
r 



/ 2m \ Q , o, , 5m, 
(1 -fc^c'r^ (1-1- 3m' ) 



+ c'c'r^ - c"cr^ 



null geodesies, A:" 



where A is an affine parameter. 



The condition for conformal flatness of the models is 
therefore W(r, v) = 0. 
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